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ABSTRACT 

We present a light curve model of the symbiotic nova PU Vul (Nova Vulpeculae 1979) that shows a 
long-lasted flat peak with no spectral indication of wind mass-loss before decline. Our quasi-evolution 
models consisting of a series of static solutions explain both the optical flat peak and ultraviolet (UV) 
light curve simultaneously. The white dwarf mass is estimated to be ~ Q.6Mq. We also provide a 
new determination of the reddening, E{B — V)= 0.43 ± 0.05, from UV spectral analysis. Theoretical 

light curve fitting of UV 1455 A provides the distance oi d = 3.8 ± 0.7 kpc. 

Subject headings: binaries: symbiotic — nova, cataclysmic variables — stars: individual (PU Vul) — 
ultraviolet: stars — white dwarfs 



1. INTRODUCTION 

PU Vul was independently discovered by Y. Kuwano 
()Kozai Ill97"9al ) and M. Honda ()Kozai Ill979bl ) as a nova- 
like object. Subsequent observations revealed that PU 
Vul did not behave like a classical nova in its photomet- 
ric and spectroscopic properties. After the rise to the 
optical maximum, PU Vul maintained an almost stable 
maximum oiV — 8.6 from 1979 to 1987 except for a deep 
minimum in 1980. It started to fade in 1988 very slowly 
toward the prediscov ery magnitude of B = 14.5 — 16.6 
(|Liller fc Liller 1 119791 ) and the photographic magnitude 
TTipg = 15.0 — 16.5 (|Yamashita et al. 1 11982D . Such a 
long-lasted flat peak and a very slow evolution of the 
light curve made this object quite different from ordi- 
nary novae. PU Vul was recognized to be a binary system 
consisting of an M giant and an outbursting co mponent 
(jBelvakina et al. lll982aHFriediung et al. II1984D and also 
that the outburst is a thermonuclear run away event on a 
white dwarf (WD) in a symbiotic binary (jKenvon II 19861 ). 
There were debates on the origin of the deep m inimum 
occurred in 1980 (see discussion in lKenvoinil986D . After 
the second eclipse in 1994, it has been clear that PU Vul 
is an eclipsing binary with the orbital period of 13.4 yr. 
The observational properties are summarized in Table 
[U that shows the discovery date, nova speed class, re- 
markable property of the light curve, indication of dust 
formation, observational evidence of eclipse, orbital pe- 
riod, peak magnitudes of my, information of UV 1455 

|mariko@educ.cc.keio.ac.jp| 



A continuum-band light-curve obtained from fitting of 
Model 1 (Section 14. 1|) . i.e., peak value and full width 
of half maximum (FWHM) , extinction in literature and 
our estimates, distance to PU Vul in literature and our 
estimates. 

Spectral development of P U Vul was exte n sively 
studied by various autho rs (jYamashita et al. ' 11982 



1983 ; lliiima fc OrtolanQ 



198S : iGochermann 



Kanamitsu et al 



Klein et al. 



^ 



. .1991! 
1991bl : 



1984; 



Belvakina et al, 
Kana mitsu I 11991a : 
iTamura et al. I 11992 : 



19941). The spectra mirnicked those of an F 



supergiant (jYamashita et al. 1 119821 : iKanamitsu I ll991aD 
in th e early phase and changed to Ao (jBelvakina et al. I 
[19891: IVoeel fc Nussbaumer I Il99l from 1983 to 1986 
as the excitation temp erat ure gradually in c reased 
([Kanamitsu et al. I ll991bD . lYamashita et al. I (|1982[ ) 
commented that the eruption must have been quite soft 
because no evidence of shell ejection, both in emission 
lines and shell absorption lines, was detected. The 
optical spectrum was strongly absorption-dominated 
until 1985 but changed t o a distinct nebular spectrum in 
the sec ond half of 1987 (jlijima II1989I : IKanamitsu et al. 1 
Il991b( ). In 1990 the star had shown rich emission lines 
in the optical and UV spectrum, which are typical in 
the nebular phase and associated with an ex t ended 
atmosphere of a WD (jVogel fc Nussbaumer I 119921 : 
IKanamits u et al. Ill991bl : iTomov et al. 1991 ). 

It is very interesting that there is no indication of 
strong winds in PU Vul in contrast to many other 
classical novae. Instead, optically thin mass-ejection 
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TABLE 1 

Observational Properties of PU Vul 



subject 



data 



comments 



discovery date 


o.oz April iy/y 


U 1 ^^za^ l^^^^aj 


nova speed class 


very slow 




light curve 


flat peak 




dust 


no 




eclipse 


yes 




orbital period 


.. 13.4 


yr 


my of flat peak 


.. 8.6 


mag 


peak of UV 1455 A flux 


.. 6.0 X 10-1=^ 


erg cm~^ s~^A~^ this work: Section|4] 


FWHM of UV 1455 A 


.. 4.4 


yr this work: Section |4] 


E{B - V) 


.. 0.29-0.5 


see references^ 


E{B - V) 


.. 0.43 ±0.05 


this work: Sectionl2.1l 


distance 


.. 1.6-7 


kpc see references ^ 


distance 


.. 3.8 ±0.7 


kpc this work: Section 14.41 


" E{B -V) = 0.4 fBclyakina 


ct al. ' 1982bT, 0.49 fFriodiune et al. II1984I). 0.4-0.5 IKenvon II1986D. 


0.50 IGochermann 1991), 0.4 


(Vog;cl &c Nussbaumcr 


1992), 0.29 ILuna & Costa II2005D 



°_d = 5-7 kpc IIBelvaki na ct al. 1982b), 5.3 kp c IBelvakina et al. ifTgsl), < 5.6 kpc IGochermann I 
[T99IT) , 1.6-2.0 kpc fVoEcl & Nussbau mcr II1992I) , 2.5 kpc l IHoard et al~lll996l) 



from W D photosphere was sugge s ted from P Cygni hne- 
profiles (IBelvakina et al. 1 Il989t fVogel fc Nussbaumer I 
[l99l ISion et al. 1119931 INussbaumer fc Vogellll996D or 
triple structure of IR emission lines (jBensammar et aLl 
Il991f) . The line width corresponds to 1100 -1200 km s~^ 
average full widths at zero intensity (jTomov et al I 



1991'), ~ 2600 km s ^ in Balmer emission wings fliiimaj 
1989} , and 550-600 km s^^ in UV spectra (Sion ct al.U 
19931 ). Mass ejection was also suggested from X-ray emis- 



sion detected with ROSAT on 10-12 November 1992 UT 
and interpreted as thermal bremsstr ahlung with a tem- 
perature of 0.22 keV (2.6 x 10^ K) (IHoard et al. I [1991 
IMiirset et al. |[T997[ ). These authors suggested a colli- 
sional origin of the X-ray between a high-density, low- 
velocity cool wind from the M giant and a low-density, 
high-velocity hot wind from the WD in the context of 
common properties of symbiotic novae. 

To summarize, the outburst of PU Vul was very quiet 
in the first ten years, and optically thin mass-ejection 
arises from 1988-1990 as the nova entered a coronal 
phase. These spectral features as well as the long-lasted 
flat peak make PU Vul quite different from many other 
classical novae in which optical magnitude decays quickly 
from its peak and spectrum indicates strong optically 
thick winds. This paper aims to model such a quite dif- 
ferent evolution of PU Vul and to understand the cause 
o f such properties. 

iKato fc Hachisu I (|2009[ ) reexamined the conditions of 
occurrence of optically thick winds and found that opti- 
cally thick winds occur in a limited range of the envelope 
(ignition) mass. For a relatively large envelope mass, op- 
tically thick winds are suppressed in a way that a large 
density-inversion layer appears and the gas-pressure gra- 
dient balances with the radiation-pressure gradient, the 
driving force of the winds. In massive WDs (> 0.7 Mq), 
optically thick winds always occur, because the igni- 
tion mass of the wind-suppressed solutions are too mas- 
sive to be realized in the actual novae. In less mas- 
sive WDs (< 0.5 M0), on the other hand, no winds 
are accelerated because the radiation-pressure gradient 
is too weak to drive the winds. In between them, i.e., 
0.5 Mq < MwD < 0.7 Mq, both types of solutions 
(wind and wind-suppressed) can be realized depending 



on the initial envelope mass. For a less massive enve- 
lope, optically thick winds are accelerated and a shell 
flash develops as a normal nova with strong winds. On 
the other hand, if the initial envelope mass is relatively 
large, optically thick winds are suppressed and a nova 
e volves without winds. 

IKato fc Hachisu I (|2009D also presented an idea that 
such wind-suppressed (no optically thick wind) evolution 
will be realized in a very slow nova that shows a long- 
lasted fiat optical peak. If the optically thick wind oc- 
curs, as in many classical novae, the strong winds carry 
out most of the envelope matter in a short timescale, 
and the optical brightness quickly decays. Thus, the 
light curve has a sharp optical peak. On the other hand, 
in the wind-suppressed evolutions, the brightness decays 
very slowly, because the evolution timescale is deter- 
mined only by hydrogen nuclear burning. Therefore, the 
nova stays at a low surface temperature for a long time, 
which makes a long-lasted flat optical peak. PU Vul is 
the first example of this new type of evolution. 

In Section [2] we review observational results based on 
the lUE spectra. In Section |3| we briefly introduce our 
method and assumptions of the theoretical model. Sec- 
tion m shows how to estimate the WD mass from light 
curve fittings. Discussion and conclusions follow in Sec- 
tions [5] and [6l respectively. 

2. UV OBSERVATIONS 

PU Vul had been monitored by lUE from Febru- 
ary 1979 to September 1983 and from October 1987 to 
September 1996 at both low and high resolutions. A 
gallery o f UV SWP spectra fr o m 19 92 to 1995 can be 
found in 'Nuss baumer fc Vogel I (11996). 

In the following we revisit the problem of the color ex- 
cess E{B — V) of PU Vul, and describe the long term 
evolution of the emission lines. The evolution of the 
UV continuum will be described in Section 14.11 The 
ultraviolet spectra were retrieved from the lUE archive 
through the INES {lUE Newly Extracted Spectra) sys- 
terrQ, which also provides full details of the observations. 
The use of lUE INES data is particularly important for 
the determination of reddening correction because of the 

http://sdc.laeff.inta.es/ines/ 
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implementation of upgraded spectral extraction and flux 
calibration procedures compared to previously published 
UV spectra. 

2.1. Reddening Correction 

Table [T] shows that the color excess E{B — V) toward 
PU Vul in the literature lies in the range from 0.3 to 0.5. 
Given the large spread of these determinations, we have 
directly determined E{B — V) from the strength of the 
2200 A feature seen in the UV sp ectra of PU Vu l. 

The Galactic extinction curve (jSeaton I I1979D shows 
a pronounced broad maximum around 2175 A due to 
dust absorption. Since it takes the same value X{X) = 
A{X)/E{B - T/) « 8 at A = 1512, 1878, and 2386 A, the 
slope of the straight line passing through the continuum 
points at these wavelengths is insensitive to E{B — V) 
in a (A,logi^(A)) plot. This circumstance can be used 
to get a reliable estimate of E{B — V) as that in which 
the stellar continuum becomes closely linear in the 1512- 
2386 A region, and passes through the continuum points 
at the above wavelengths. From 9 pairs of short and long 
wavelength lUE spectra taken from JD 2,448,217 to JD 
2,450,342, i.e. during the nebular phase, we have in this 
way found E{B - V) = 0.43 ± 0.05. Examples of lUE 
spectra of PU Vul corrected with E{B — V)— 0.43 are 
reported in Figure [TJ 

2.2. Evolution of the UV Continuum 

We have measured the mean flux in two narrow 
bands 20 A wide centered at 1455 A and 2855 A, 
selected to provide a fairly good representation of 
the UV co ntinuum because little afi^ected by emis- 
sion lines (jCassatella et al. I I2002D . Figure [5] shows 
the time evolution of the F(1455 A) and i^(2885 A) 
fluxes and of the UV color index C(1455 - 2885) = 
-2.5 log[F(1455 A) /F(2885 A)] . The measurements were 
made on well exposed low resolution large aperture spec- 
tra. Figure[2]reports also, for comparison, the visual light 
curve obtained from the Fine Error Sensor (FES) counts, 
T^FES, on board lUE, once corrected fo r the time depen- 
dent sensitivity degradation (p. 35 in iCassatella et al. I 
1200 4V Note the pronounced UV maximum around JD 
2,447,400 followed by a progressive decay at the same 
time as the UV spectrum becomes harder, as indicated 
by the decrease of the UV color index. The delay 
of the UV maximum with respect to the visual max- 
imum is common both to novae and symbiotic star s 
(|Fernandez- Castro et 3171119951: ICassatella et al. II2002D . 

2.3. Evolution of the UV Emission Lines 

Figure |3] reports, as a function of time, the observed 
fluxes in the emission lines of C III] 1909 A, C IV 1550 A, 
N IV] 1487 A, He II 1640 A, and N V 1240 A, which cover 
a wide range of ionization conditions (the corresponding 
ionization energies Xion are 24.4, 47.9, 47.4, 54.4 and 77.5 
eV, respectively) . The flux measurements were obtained 
from well exposed lUE low resolution spectra. The fig- 
ure shows clearly that the eclipse around JD 2,449,550 
endures longer for high than for low ionization lines. This 
suggests that the high ionization lines are formed in the 
unseen side of the cool giant's wind, so confirming the 
results bv .Nussbaumer fc Vogel ■ (|1996i) (see their Figure 
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Fig. 1. — lUE spectra of PU Vul obtained at different dates. 
Fluxes are in units of erg cm~^ s"'^ A~^. Tfie spectra have been 
corrected for reddening using E(B — V) = 0.43. The vertical dotted 
lines represent the wavelengths AA 1512, 1878 and 2386 A at which 
the extinction law takes the same value. With the adopted value 
of reddening, the stellar continuum underlying the many emission 
lines is well represented by a straight line all over the full spectral 
range. Saturated data points in the emission lines are labeled with 
pluses. 

5), who found that the UV highest excitation lines of 
He II A1640, N V A1240, and N IV A1718 disappeared 
during the second eclipse. 

3. MODEL 
3.1. Evolution of Nova Outburst 

Nova is a thermo n uclear runaway event on a WD 

fariai et aTl [19801: iPrialnik I [19861 iStarrfield et al I 
r2, 19881 IJose et al. II1998I: IPrialnik fc Kovetz II1995D . 
After unstable hydrogen nuclear burning triggers a nova 
outburst, the envelope on the WD greatly expands to a 
giant size. After it reaches the optical peak, the envelope 
expansion settles down into a steady-state. The opti- 
cal magnitude decreases as the envelope mass decreases 
and the photospheric temperature rises with time. In 
less massive WDs the optically thick wind does not oc- 
cur as described in Section [T] and, therefore, the decay 
phase of nova e can be fol l owed by a quasi-hydrost atic se- 
quence (e.g., [ibeniTgsl IKato fc HachisuT[T99l [2009). 
We solved the equations of hydrostatic balance, radiative 
diffusion, and conservation of energy from the bottom 
of the hydrogen-rich envelope through the photosphere. 
The bottom radius is assumed to be the Chandrasekhar 



4 



Kato et al. 



2i 



15 



-5 



Flux 1455 A 




Flux 2885 A 




C{1 455-2885) 



FES Visual Flux 




180 



4000 



5000 



6000 
JD - 



7000 8000 
2440000 



9000 10000 



Fig. 2. — Evolution of the continuum fiuxes at 1455 A and 2885 
A, of the ultraviolet color index C(1455 — 2885), and of the Vfes 
visual flux of PU Vul (see Section I2.2| l. Visual and UV fluxes 
are both in units of 10~^^ erg cra~^ A~^, not corrected for 
reddening. Only the color index ha s be en corrected for reddening 
using E{B — V) = 0.43 (see Section [2ll . The vertical dotted lines 
indicate the period of the first (JD 2,444,380 - 2,444,710) and the 
second (JD 2,449,275 - 2,449,608) ecUpses. 



radius. The evolution of novae is followed by connecting 
these solutions along the envelope mass-decreasing se- 
quence. The time evolution is calculated from the mass 
decreasing rate which is the summation of the two rates, 
hydrogen nuclear burning and optically-thin wind mass - 
loss. We used OPAL opacities (jlglesias fc Rogers"1ll996[ ). 
The method and numerical techniques are essentially the 
same as those in [Kato & Hachisu (1994). Convective 
energy transport is calculated using the mixing length 
theory with the mixing-length para meter a = 1.5 (see 
Figure 11 in lKato fc Hachisul 120091 for the dependence 
of model light curve on the mixing length parameter a). 

In the rising phase, the envelope does not yet settle 
down to a thermal equilibrium, i.e., the nuclear energy 
generation is larger than the radiative loss. We have ap- 
proximated such a stage by a sequence of static solutions 
of constant mass without thermal equilibrium. These so- 
lutions may not approximate well the rising phase, but 
are enough to our purpose, because the rising phase plays 
no important role in the determination of physical values 
such as the WD mass and distance. 

3.2. Wind Mass-Loss Rate 
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Fig. 3. — Observed fluxes of some prominent UV emission lines 
of PU Vul as a function of time. Fluxes are in units of 10~^^ erg 
cm~-^ s~^, not corrected for reddening. 



In the later phase of the outburst (after 1986), 
optically-thin wind may arise because many emission 
lines had appeared in spectra. We cannot calculate such 
wind mass-loss, however, because radiative transfer in 
the optically-thin region is not included in our model 
calculation. Therefore, we take the wind mass-loss rate 
as a model parameter. This optically-thin wind does not 
affect the envelope structure below the photosphere, but 
speeds up the nova evolution because mass decreasing 
rate of the envelope is accelerated by this wind mass- 
loss. 

In line driven winds, the mass-loss rate M is limited 
by photon momentum; the wind cannot get momentum 
much exceeding the momentum of photon flux. This is 
in contrast to optically thick winds in which the mass- 
loss rate co uld be much larger than the momentum of 
photon flux ()Kato fc Ibenlll992D . Therefore, we assume 
a condition, that the momentum of wind is smaller than 
that of photon flux, of 



MWwind < 



(1) 



where c is the speed of light. We get an upper limit of 
the wind mass-loss rate as 
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M < 1.6x10"^ ( 



L 



ph 



^wind 



6 X lO^^erg s-i'"'200 km s 



Here we take Lph = 6 x 10^^ erg s^"'^ from a bolometric 
luminosity of a 0.6 M© WD (which will be shown later 
as Model 2 in Table [J) , and assume a relatively small 
v alue of fwinri = 200 km for a safe upper limit. 

iBensammar et aL~l ()1991t ) estimated the mass-loss rate 
of PU Vul to be about 4 x 10"^ Mq yv'^ from emission 
measure on 30 April and 1 May 1988 with a wind ve- 
locity of 70 km s~^ and a hot star radius of 3.9 x 10^^ 
cm (56 Rq). iSion et al~l (|1993[ ) set an upper limit of the 
mass-loss rate, 1 x lO'^ Mq yr'^. iSkopal I (|2006D esti- 
mated the mass-loss rate from Ha line luminosity of 144 
Lq on 28 September 1988, assuming optically-thin fully- 
ionized winds. With a terminal velocity 2,100 km s^^ 
of a given velocity profile, Skopal obtained the mass-loss 
rate to be Af = 5 x 1O~^M0 yr~^. These values, however, 
may be a little bit larger than our momentum condition 

In our model, we assume no wind mass-loss during the 
optical flat peak because no emission lines are observed 
or they are very weak. We assume that optically thin- 
wind begins when the photospheric temperature rises to 
logTph (K) ~ 4.0. This wind seems to have much weak- 
ened sometime around 1999-2001, because, after that, 
the optical light curve changed its shape and seems to 
decay as t~^, where t is the time after the decay started. 
In classical novae, the decay is often observed in the 
later phase of the outburst, and is interpreted as emission 
from homologously expanding optically thin plasma with 
a constant m ass, which indicates that the mass supply 
had stopped pachisu k Kato 1120061 120101) . In the case 
of PU Vul, the wind had not stop ped entirely since P Cyg 
profiles are stih observed in 2004 ()Yoo Il2007[) . Therefore, 
we assume that the optically-thin wind begins at log Tph 
(K) = 4.0 and continues until log Tph (K) — 5.05 at a 
rate shown later, and after that the wind mass-loss rate 
dropped to M = 1.0 x lO^'^Af© yr"^ 

3.3. Multiwavelength Light Curves 

After the maximum expansion of the photosphere, the 
photospheric radius (Rph) gradually decreases keeping 
the total luminosity (iph) almost constant. The photo- 
spheric temperature (Tph) increases with time because of 
Lph = 47ri?pj^crrp[j. The maximum emission shifts from 
optical to supersoft X-ray through ultraviolet (UV). This 
causes the luminosity decrease in optical and increase in 
UV and finally increase of supersoft X-ray. We assume 
that photons are emitted at the photosphere as a black- 
body with a photospheric temperature Tph. The light 
curve of optical (V) and UV 1455A fluxes are estimated 
from the blackbody emission. 

3.4. Chemical Composition 

'Bclvaki na et al. 1 (|1989l ) obtained chemical composi- 
tion of the atmosphere of the hot component of PU Vul, 
in which iron is depleted by a factor of 0.3-0.5 against the 
sun. The n umber ratio of He/H is estimated to be 0.31 
(jAndrillat fc Houziaux ..1994) and 0.146 (Luna fc Costa I 



[20051) from emission line ratios, which shows helium over- 
abundance than the solar value (He/H ^ 0.08). An- 
-)~^ Mq yr~^. other suggestion comes f rom the location in the Galaxy. 

IBelvakina et al. 1 (jl9821i ) suggested that PU Vul does not 
(2) belong to the planar component of the Galaxy because 
the star is off the galactic plane by 0.7-1.0 kpc, based on 
their derived distance of 5-7 kpc. This value, however, 
reduced to 0.5 kpc if we adopt the distance of ~ 3.8 kpc 
as we will obtain later. 

With the information above we assume the mass frac- 
tion of hydrogen, helium and heavy elements of the enve- 
lope to be (X, Y, Z)=(0.5, 0.494, 0.006). For compari- 
son, we further assume additional sets of composition i.e., 
(0.5, 0.49, 0.01), (0.7, 0.28, 0.02) and (0.7, 0.29, 0.01). 
We simply assumed that the chemical composition of the 
envelope is uniform and constant with time. 



4. LIGHT CURVE FITTING 

4.1. UV Light Curve Fitting and the WD mass 

The UV 1455 A flux provides a good representation 
of the continuum level in novae, bec ause it coincides 
with a local minimum of line o pacity (ICassatella et al. 1 
^00^. In previous papers CHac hisu fc Kato I 120061 



Hachisu et aLH 2008: Kato et al. 200^ we have shown 
that the UV 1455 A continuum light curve is very sen- 
sitive to model parameters, especially the WD mass. 

Figure 0] shows time evolution of the UV 1455 A con- 
tinuum fiux as well as our theoretical light curves. Figure 
|4^ shows dependence on the WD mass for a given set of 
the wind mass-loss rate of 1 x lO^^M© yr~ (arbitrar- 
ily chosen but not too small compared with the nuclear 
burning rate of 3 x 1O~^M0 yr~^) and chemical compo- 
sition of X = 0.5 and Z = 0.006. We see that the UV 
peak is narrower in more massive WDs, because nova 
evolves faster owing to a less massive envelope and a 
high nuclear burning rate. Figure |3)d depicts four light 
curves with different wind mass-loss rates for a given WD 
mass and chemical composition. As the assumed wind 
mass-loss rate is comparable to the nuclear burning rate 
(3 X 10~^Mq yr~^), the evolution speed of nova is sen- 
sitive to the mass-loss rate. Figure Hj; shows five light 
curves with different sets of chemical composition. Nova 
evolves faster for smaller X, because of less nuclear fuel, 
and also faster for larger Z because of a smaller envelope 
mass. 

In this way we can choose a WD mass with reason- 
able agreement with the UV data for a given parameter 
set of the wind mass-loss rate and composition, {X, Y, 
Z)=(0.5, 0.494, 0.006). Table [D shows two such mod- 
els, one is a 0.57 Mq WD with the wind mass-loss rate 
of 5 X IO-'^Mq yr-i (Model I) and the other is a 0.6 
Mq WD with 3 X 10-'^ Mq yr'^ (Model 2). These WD 
masses, 0.57 Mq and 0.6 Mq, are not much different, 
because they are in the middle of the permitted range of 
the WD mass, 0.53-0.65 Mq, as explained below to avoid 
extremely small and large mass-loss rates. The 0.57 and 
0.6 Mq models correspond respectively to the left and 
upper sides to the "wind region" (a triangle region in 
Figure 10 of lKato fc Hachisiol 120091 ) for the correspond- 
ing composition. 

During the outburst, the photospheric temperature 
gradually rises and the photospheric radius decreases, 
with an almost constant photospheric luminosity. Figure 
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TABLE 2 

Summary of Our Models 



subject 


model 1 


model 2 


units 


X 


. 0.5 


0.5 




Y 


. 0.494 


0.494 




Z 


0.006 


0.006 




WD mass 


. 0.57 


0.6 


Mq 


distance from UV fit ^ 


. 3.7 


3.9 


kpc 




. -5.49 


-5.58 


mag 


r b 

^pcak 


. 5.2 


5.7 


lO'^^erg s~l 


H-burning rate 


. 2.7 


3.0 


lO-^A/oyr-l 


maximum radius 


. 62 


63 


Rq 


initial envelope mass 


. 5.8 


4.6 


IO-^Mq 


assumed wind mass-loss rate (T < 5)° 


. 5.0 


3.0 


10"^ Mq yr-i 


assumed wind mass-loss rate (T > 5)' 


. 1.0 


1.0 


10"^ Mq yr~i 


mass lost by the wind (T < 5)° 


. 0.60 


0.34 


IO-^Mq 


mass lost by the wind (T > 5)* 


. 0.61 


0.44 


lO-^AfQ 



with E{B -V) = 0.43 
^ values at log Tph (K)=3.9 
values at log Tph (K)=4.5 



the mass at the rising phase. 
^ optically-thin wind from log Tph (K) — 4 to 5.05. 
^ optieally-thin wind at log Tph (K) > 5.05 

[5] shows the developments of the temperature and radius 
of Model 1 and Model 2 as well as observation al estimates 
taken from literature. iVogel fc Nussbaumer I (j.1992) sug- 
gested that the temperature was as low as 6,000 K in 

1979 as an effective temperature of F-supergiant and 
estimated the temperature to be 40,000 K < T < 
50,000 K at 1991 from nebular spectra. They also ob- 
tained t he radius of ^ 5 Rp, from the ana l ysis o f 1980 
eclipse. iKenyonI (I1986D . iBelvakina et al. I (|1989f ). and 
IGochermann I (jl99ll ) estimated the temperature from 
spectral type. These values are plotted in Figure [SJ in 
which our model temperature well reproduces observa- 
tional estimates including the epoch around the UV1455 
A light curve peak, i.e., JD 2,447,000 to 2,448,200. The 
model radius is also consistent with observations consid- 
ering large ambiguity in the estimating methods. 

Table [3] shows the upper and lower limits of the WD 
mass that reproduces reasonable fitting to the UV light 
curve. For a given chemical composition, the maximum 
WD mass is obtained with no wind mass-loss, because 
more massive WDs produce much narrower UV curves. 
A more massive WD than the maximum value in Table 
[3] produces a UV light curve too narrow to fit the ob- 
servation. On the other hand, for less massive WDs, we 
need to assume larger wind mass-loss rates because of 
their slower evolutions. The minimum WD mass may be 
obtained for a largest wind mass-loss rate that we adopt 
1 X 10~^Mo yr~^ from Equation ([2]). With plausible val- 
ues of the optically thin wind-niass-loss rate, which may 
be a few to several x10~^Mq yr~^, we may conclude 
that the WD is about 0.6Mq. 

4.2. Optical Light Curve 

Figure [6] shows optical light curves of Model 1 and 
Model 2, of which characteristic values are summarized 
in Table [21 These two models are selected from fitting 
with the UV light curves, but also well reproduce the 
optical light curve in the flat maximum as well as the 
following decline until 1989, except the first eclipse in 

1980 which is not taken into account in our model. 
After 1989 the theoretical light curve largely deviates 



from observed optical magnitudes. In this stage, spectra 
are e mission-line dominated and the continuum is very 
weak (|Yooll2007[ ). These emission lines comes from op- 
tically thin plasma outside the photosphere which is not 
included in our model. Therefore, our theoretical models 
give much lower magnitudes than that of observational 
data. 

In our theoretical models, the flat peak corresponds 
to the era of low photospheric temperature (7,000-9,000 
K) as shown in Figure [3 The temperature gradually 
rises with time and reaches 10,000 K which is indi- 
cated by the cross in Figure [SI where we assume the 
optically-thin mass-loss begins. After that the nova en- 
tered t he coronal phase and many emission lines ap- 
peare d (jKanamitsu et al.1ll99lbl : INussbaumer fc Vogel I 
11996( 1. Our model temperature is consistent with these 
observational properties. 

4.3. Internal Structure of the Envelope at the Flat Peak 

Figure [7] shows internal structures of the envelopes of 
Model 1 and Model 2 with the photospheric temperature 
log Tph (K) ^ 3.9. Here, the local Eddington luminosity 
is defined as 

_ 47rcGAfwD 

^Edd = , (o) 

K 

where k is the opacity in which we use the OPAL opac- 
ity. Since the opacity k is a function of temperature and 
density, the Eddington luminosity is also a local vari- 
able. This Eddington luminosity has the deepest local 
minimum at logi? (cm) — 10.4 -11.2 that corresponds 
to the Fe peak of OPAL opacity. There appears a large 
density-inversion layer at logr (cm) ~ 10.2 — 11.2 corre- 
sponding to the super Eddington region (l^Edd < L^)- 
This density-inversion arises in order to keep hydro- 
static balance in the super-Eddington region. Such a 
structure is very different from ordinary nova wind so- 
lutio ns, in which the density mon otonically decreases as 
|Kato fc Hachisu Ill994f . l2009f ). but similar to that of 
red giants. 
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TABLE 3 

Upper and Lower Values of the WD Mass 







composition 




minimum mass^ 


maximum mass'' 


X 


= 0.5, 


Y = 0.294, Z -- 


= 0.006 . 


. 0.53 Mq 


0.65 Mq 


X 


= 0.5, 


Y = 0.29, Z = 


0.01 


. 0.5 Mq 


0.62 Mq 


X 


= 0.7, 


Y = 0.28, Z = 


0.02 


. 0.5 Mq 


0.67 Mq 


X 


= 0.7, 


Y = 0.29, Z = 


0.01 


. 0.52 Mq 


0.72 Mq 



in case of a very large mass-loss rate of 1 x 10 ^ Mq yr ^ 



extreme case of no wind mass-loss 
3 



— 1 ' r 

UV1455A a 

1.E-7 Mq yr 
X=0.5 
Z=0.006 




6000 



8000 

JD 2,448,000+ (day) 



Fig. 4.— UV 1455 A light curves. The UV 1455 A data are 
denoted by red open circles (outside the eclipse) and red filled 
circles (during the eclipse). Theoretical UV light curves are shown 
in curves for an assumed distance of 1 kpc (with right-side axis), a) 
Comparison among different WD masses for the same wind mass- 
loss rate (1 X 10~ Mq yr~') and the same chemical composition 
of X = 0.5 and Z = 0.006. The WD mass is indicated beside each 
line, b) Comparison among different mass- loss rates for the 0.57 
Mq WD with composition of X = 0.5 and Z = 0.006. The wind 
mass- loss rate (in units of Mq yr~') is attached to each curve, c) 
Comparison among different chemical compositions for the same 
mass-loss rate of 5 X IO^^Mq yr"-*^ and WD mass of 0.6 Mq. 
The chemical composition is (X, Z)= (0.7,0.004) for a three-dot- 
dashed curve, (0.5,0.006) for a thick solid curve, (0.7,0.01) for a 
thin solid curve, (0.5,0.01) for a dash-dotted curve, and (0.7,0.02) 
for a dashed curve. 

These two solutions are representative of those in the 
long-lasted flat-peak in the optical light curve (Figure H]), 
i.e., the envelope is extended to 50-60 Rq and the tem- 
perature is as low as logTph (K) ~ 3.9-4.0. The convec- 
tion, which is dominant in energy transport in the rising 



1980 



1985 



1990 




4000 



6000 

JD 2440000+ (day) 



8000 



Fig. 5. — Development of the temperature and radius of the 
hot component (WD photosphere) in PU Vul. Curves denote the 
blackbody temperature and the photospheric radius of our Model 1 
(Solid line) and Model 2 (Dashed line). The effective temperature 
(open symbols) and radius (filled symbols) estimated by various au- 
thors are also sh own. Circle s: Vogcl & Nussbaumer (1992). Tri- 
angles: .Bel vakina et al. I I ll989i ). Diamonds: iGochermann i (|1991| ). 
Square: IKenvon I 1 1198^ 



phase of nova outbursts (|Prialnik II1986D . has retreated 
and is ineffective in and after the flat-peak. The con- 
vection occurs where the opacity decreases outward, i.e., 
the local Eddington luminosity increases outward. The 
largest convective region is at logr (cm) = 10.5-11.4, 
corresponding to the super-Eddington region (see Figure 
[7|). In all convective regions, convection is ineffective in 
energy transport due to low density, and unable to carry 
all of the energy flux. Therefore, the structure is super- 
adiabatic, i.e., entropy decreases outward. This situa- 
tion is different from the convective core of intermediate- 
mass main-sequence star or the inner convective enve- 
lope of low-mass red giant star, where the convective en- 
ergy transport is effective and the tempe rature gradient 
is very close to the adiabatic gradient (jHavashi et al. I 
1962t'lben 1965 ;). 

The Eddington luminosity with electron scattering 
opacity LEdd.ei = 47rcGM/[0.2(l-|-X)] is often used as an 
easy estimate of the WD luminosity. We note, however, 
that the photospheric luminosity is only 54 % of L^dd ci 
for Model 1 (0.57 Mq WD) and 57 % for Model 2 (0^.6 

Mq WD). 

4.4. Distance 
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... 1 : . — - • — ■ 



5000 10" 1.5x10" 



JD 2,440,000+ (day) 

Fig. 6. — Optical and UV light curves of PU Vul. Optical 
data are taken from lAU Circular (for dip: 3421, 3477, 3494, 
3589, 36 04, 3610, and 36 55), Wcnzcl (1979), Yamashita ct al. 
(1^85), 'Belvakina et al. ' (1982b), Chochol. Hric. & PapousekJ 
©81), Pureathofcr & Schncll (1982), Kolotilov (1983)7 
KrEathofcr & Schnell (1983), liiima & Ortolani (1984), 
lijima (1989), Kaiiamitsu ct al. (1991b), Klein ct al. (1994), 
Kolotilov et al. , a995i) . lYoon fc Honevcutt . 12000.) . and AAVSO 
(after JD 2,452,000). Large red open circles denote the flux of 
IVE UV 1455 A band. Calculated light curves are also shown for 
Model 1 (Thin solid line) and Model 2 (Dashed line). The scale in 
the right-hand-side axis denotes that for observational data. Scale 
for the theoretical UV flux is (bottom, top)=(0, 6.25) for Model 1 
and (0, 7.0) for Model 2 in units of 10"'' erg s"! cm'^ A"! for 
the distance of 10 pc. The cross/small dot indicate starting/end 
point of the assumed optically-thin wind mass-loss. See text for 
more details. 

The distance to PU Vul is obtained from the compar- 
ison of the 14 55 A band flux with t he corresponding 
model fluxes (jHachisu fc Kato 1 120061 : IKato fc Hachisu I 
[200I [2007I) . The flux of Model 1 is = 
2.2 xlO"*^ (d/10 pc)2 ergs cm^^ g-i ^-1 ^j^g 
peak. The corresponding observed flux is F'^^^—QSi 
xlO~^^ ergs cm~^ s~^ as in Figure [51 From these 

values we obtain the distance of d = 3.7 kpc wit h the ab- 
sorpt ion Ax = 8.3 X E{B - V) for A = 1455 A (jSeaton I 
Mm . here we use E{B ~ V) — 0.43 (obtained in 
Section I2.ip . If we use the upper and lower limit of 
E{B~V) = 0.43±0.05, we get the distance oid = 3.7tol 
kpc. The error coming from the UV flux fltting is much 
smaller and up to ±0.15 kpc. In the same way, we ob- 
tain d — 3.9^0 6 foi' Model 2. If we adopt a different set 
of chemical composition, the resultant distance is also 
changed. For the 0.6 Mq WD we obtain 3.8 kpc for 
(X, Z) = (0.5, 0.01), 3.3 kpc for (0.7, 0.02) and 3.5 
kpc for (0.7, 0.01). Taking into account that Z = 0.02 
is a bit larg er than the rec ent estimate of solar value 
{Z = 0.0128: iGrevesse | I2008). 3.3 kpc may be a smaller 
limit. Thus, we may summarize our distance estimates 
as d = 3.8 ± 0.7 kpc. The error includes ambiguity of the 
E{B — V), WD mass, and chemical composition. 

5. DISCUSSION 

As described in Section [3l a nova becomes a super- 
soft X-ray source in the later phase of the outburst. In 
PU Vul, however, the maximum temperature is not high 
enough compared with a typical classical nova because of 
the less massive WD 0.6^/0). Our theoretical model 
predicts that the outburst of PU Vul is still ongoing and 




log r (cm) 

Fig. 7. — Temperature (T), density (p), diffusive luminosity (L^), 
and the local Eddington luminosity (LEdd) of envelope solutions 
at logTph = 3.9 for Model 1 (solid curve) and Model 2 (dotted 
curve). The right edge of the curve corresponds to the photosphere 
and the left edge to the bottom of the envelope, i.e., bottom of 
the nuclear burning region. Two short horizontal lines indicate 
Eddington luminosity with electron scattering opacity, -^Edd el = 
47rcGM/[0.2(l -I- X)] for Model 1 (lower line) and Model 2 (upper 
line). 

the temperature is continuously rising. The temperature 
will finally reach the maximum temperature of ^ 3 x 10^ 
K just before hydrogen burning stops (X-ray turnoff), 
while the flux is almost constant at 2 — 3 x 10"^^ erg s~^. 
Therefore, we may expect supersoft X-ray from PU Vul 
in future, but it is very difficult to predict this epoch, 
because it strongly depends on the model parameters, 
such as the optically-thin wind mass-loss rate, and the 
chemical composition of the envelope. 

No X-ray observations of PU Vul have been made since 
November 1992. The detected X-ray flux was attributed 
to a shock origin oi^ colliding winds (jHoard et al. 1 119961 : 
IMiirset et al~ni997l) . PU Vul is a symbiotic nova in 
which the M-giant companion blows a massive cool wind 
that may preferentially distribute in the orbital plane. 
This cool wind may prevent a clear detection of super- 
soft X-rays from the hot WD. Therefore, the probability 
of detecting supersoft X-rays would be higher when the 
WD is in front of the cool wind of the giant (around 2014) 
and would be very low near the eclipse (2020-2021). 

We have performed simulations for detectability of su- 
persoft X-ray with the assumed values of E{B — V) = 
0.43 and d ~ 3.8 kpc. It will be possible to detect PU 
Vul at 0.007 cts s^^ with the EPIC-pn camera onboard 
XMM-Newton when the temperature is logT (K)^ 5.3 
and the luminosity is ^ 5 x 10'^'' erg s~^ (probably in 
2010-2020) and at about 1 cts s~^ near the X-ray turnoff 
time, i.e., logT (K)^ 5.5 and ~ 2 x 10^'^ erg s~^ (prob- 
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ably in 2060-2090). However, these estimates are based 
on our Model 1 and Model 2 only, and there is a large 
ambiguity due to uncertainty of the model parameters 
and assumptions. Especially the above expecting year 
depends strongly on the assumed wind-mass-loss rate af- 
ter the photospheric temperature of the WD envelope 
rises to logT (K)=: 5.05. If the wind-mass-loss will much 
weaken, the X-ray turnoff time is much later than the 
above estimates. 

We expect that a large part of the envelope will re- 
main on the WD after the outburst of PU Vul, because 
no optically thick wind occurs. From the initial enve- 
lope mass and the matter lost by the optically-thin wind 
as in Table [2l we estimate the envelope mass that will 
remain after the outburst. Highly depending on the as- 
sumed mass-loss rate as well as the other parameters, 
it is estimated to be about 70-90 % of the initial enve- 
lope mass for reliable model parameters. Mass-accreting 
WDs are, in general, potential candidates of progenitors 
of Type la supernovae (SNe la). However, low mass WDs 
(< O.9M0) such as in PU Vul have not been considered 
as a candidate of SNe la, because it is hard t o grow to 
Chandrasek har mass limit (e.g. iHachisu et al. .,19998^ ^ 
[20101 ). 



6. CONCLUSIONS 

Our main results are s ummarized as follows: 

1. Based on the idea of iKato fc Hachisul (|2009( ) that a 
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long-lasted flat peak of optical light curves can be repro- 
duced by a sequence of wind-suppressed static-solutions, 
we have succeeded in reproducing the long-lasted optical 
flat peak of PU Vul as well as the UV 1455 A contin- 
uum light curve. Our model is consistent with spectral 
features with no indication of strong winds in the flat 
peak of PU Vul. 

2. An analysis of the lUE spectra of PU Vul indicates 
E{B -V)= 0.43 ±0.05. 

3. We obtain a mass range of the WD between 0.55 
and 0.65 Mq by comparing our theoretical light curves 
with the UV light curve. 

4. We obtain the distance of d = 3.8 ± 0.7 kpc with 
E{B - F) = 0.43 ±0.05. 

5. We may conclude that the outburst of PU Vul is 
still on-going, and has already entered the supersoft X- 
ray phase. We encourage X-ray observations in ^ 2014, 
when the WD is in front of the red giant, that will provide 
important information on this symbiotic nova. 
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